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^ ■ Abstract 

We show how to compute transport coefficients in gauge theories by considering the expansion 
of the Kubo formulas in terms of ladder diagrams in the imaginary time formalism. All summa- 
tions over Matsubara frequencies are performed and the analytical continuation to get the retarded 
. correlators is done. As an illustration of the procedure, we present a derivation of the transport 
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"^j- , equation for the shear viscosity in the scalar theory. Assuming the Hard Thermal Loop approxima- 

te 

tion for the screening of distant collisions of the hard particles in the plasma, we derive a couple of 
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integral equations for the effective vertices which, to logarithmic accuracy, are shown to be identical 
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to the linearized Boltzmann equations previously found by Arnold, Moore and Yaffe. 
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I. INTRODUCTION 



The development of a transport theory for QCD in the regime of high temperature has 
turned out a valuable pursuit with the advent of heavy ion colliders which provide novel 
tools for the study of the properties of highly excited matter. From a purely theoretical 
point of view, the computation of transport coefficients amounts a challenge even in weakly 
coupled theories, because these quantities usually depend non-analytically on the coupling 
constant. In most previous computations ]T|, [|, a kinetic approach based in the Boltzmann 
equation has been used. It is only recently within this framework that a reliable computation 
to logarithmic accuracy in gauge theories has been reported [|J. The complete leading order 
is still unavailable except for the case of a gauge theory with a large number of fermionic 
species ||. 

However, there exists an alternative approach based on Kubo formulas for appropriate 
correlation functions, which has been largely used in the context of low energy many-body 
physics 0. For the electrical and thermal conductivities of ordinary metals and supercon- 
ductors, the computation of the current correlators requires the resummation of an infinite 
class of ladder diagrams, a task that in the complete relativistic setting of gauge theories 
usually appears as a very hard issue, partly motivating the scant of this approach. In rel- 
ativistic transport theory, the resummation of ladder diagrams has been performed for the 
scalar theory by Jeon [7[], who proved the equivalence with the Boltzmann equation, and 
since repeated a few times M, |9]. Also, a simplified ladder summation has been performed 
for the computation of the leading-log order of the color conductivity |IIJ. Recently, Arnold, 



Moore and Yaffe [11], [T2| have performed a ladder summation in order to account for the 
effect of multiple scattering in the process of photon production from a QCD plasma. 

The purpose of this paper is twofold. First, we wish to explicitly show how to perform the 
summation of a restricted set of ladder diagrams in gauge theories within the imaginary time 
formalism of thermal field theory. In contrast to the work of Jeon |7]] , we will not use the series 
of cut ladder diagrams. Rather, we closely follow a treatment due to Holstein |l3j who, a 
long time ago, performed a ladder summation in order to compute the transport properties 
of the low-energy electron-phonon gas. In this approach, one first identifies the required 
analytic continuation of the effective vertex function entering in the current correlator, and 
then writes the integral equation for this vertex, summing all ladders. As shown below, this 
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approach exactly reproduces the correct transport equation for the shear viscosity in the 
scalar theory. 

On the other hand, we will try to derive the logarithmic accuracy of the transport coeffi- 
cients by only considering the role played by the soft degrees of freedom which are exchanged 
in the collisions between the plasma constituents. This requires the use of the Hard Ther- 
mal Loop (HTL) approximation |TJJ for the internal lines associated with the rungs of the 
diagrams, and the introduction of an arbitrary momentum scale q c separating the hard and 
soft ranges of the momentum transfer [^q| . An important step towards the complete com- 
putation of the hard contribution was already recently made by the authors of Ref . @] , who 
calculated the infrared logarithmic divergences of the collision terms of linearized Boltzmann 
equation written in terms of unscreened interactions. 

Our main results are a couple of integral equations for the effective vertices, encoding 
the effects of distant collisions in the plasma coming from soft momentum transfer q < q c . 
Although these equations are necessarily incomplete, they reproduce the required logarithmic 
dependence on q c which makes possible the eventual cancellation of the arbitrary scale q c in 
the final result. Hence, they reproduce the known results M for the transport coefficients 
to logarithmic accuracy. 

The plan of this paper is as follows. In Sec. II, we review some standard material on the 
imaginary time formalism of thermal field theory and Kubo formulas. Here, we include a 
useful summation formula over Matsubara frequencies and the procedure of analytic con- 



tinuation. In Sec. |H we show how to derive the shear viscosity of A0 4 theory. Section [IV 
deals with the simplifications which appear by summing the ladders when only the effect 
of soft momentum transfer is considered. Then, we derive the corresponding contribution 
to the transport equations for the electrical conductivity and the shear viscosity. Section [V| 
presents a brief derivation based on sum rules of the logarithmic terms in the transport co- 



efficients, and Sec. [VI] closing the paper contains a summary and prospects. There are short 
appendixes with some details about spectral densities, sum rules and the relevant thermal 
widths to be included in the propagators of hard particles. 
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II. BASIC FORMALISM 



A. Single particle spectral densities 

The basic element of a diagram in the Imaginary-Time Formalism is the Matsubara 
propagator depending on the purely imaginary frequencies iuj n = inn/ (3 (with n even for 
bosons and odd for fermions), 

a / • v [°° dq° p(Q) 

A(iu n ,q) = / 5 — , 1 

J-oo 2n g u — i<jj n 

where the real quantity p(Q) (with Q^ 1 = (g°, q), q° real) is the single particle spectral density. 
The analytical continuation iu n — > z defines a function A(z, q) of a complex variable z which 
is analytical off the real axis and the discontinuity through the branch cut along Imz = is 
proportional to p{Q). The different Green's functions for real frequency can be constructed 
from the spectral density. For instance, the bosonic Wightman functions A >,< (Q) are given 
by 

A>(Q) = (l + n b (q )) p(Q), A<(Q) = n b (q°)p(Q), (2) 

where rib(q°) = l/(e f3q ° — 1) is the bosonic occupation number. The retarded and advanced 
Green's functions, which will play an important role in our discussion, are 

A ret (Q) = A(g° + *0 + ,q), A adv (Q) = A(g°-*0 + ,q), (3) 

and p(Q) = 2ImA rct (Q). 

For a free particle, the spectral density is given by a superposition of delta functions, 
with support on the mass shell, p% = e 2 p . If the interactions are weak, a delta function can be 
replaced by a Lorentzian with a small width T p , related to the imaginary part of the energy 
at the poles of the retarded propagator, p° = ±e p — iT p , 

^^" Ti^ (4) 

which gives rise to an analytical propagator 

_ dp 2F p 1_ -1 

{z, pj - J ^ 2 ^ (po -p £p ) 2 + F 2 po _ ^ - z ^ £p + iY p S ign(lm z) ' U 

Thus, for this case, the function A(z, p) has a branch cut at real axis in the complex z-plane 
and it has no poles. 



At high temperature, the particles entering in the collision processes occurring in the 
plasma are mostly particles propagating nearly on-shell with hard momentum, P ~ T. 
Their spectral densities may be approximated by a combination of two Lorentzians 

Pb(P) = p {(p°-p) 2 + ti ~ (p°+pf + ri) ' (6) 

where h±(p) = (7 =p 7 • p)/2, and the thermal widths T p and 7 P are the imaginary parts of 
the transverse piece of the on-shell gluon self-energy and the quark, respectively 

r p = ~1mIE*(jf i = p,p), (8) 

7p = -ltr(^Im£ rc V=p,p)). (9) 

The shift in the real part of the energy can be ignored since it is perturbatively small when 
the energy is 0(T). 

In gauge theories, the imaginary part of the thermal self-energies receives contributions 
from various scattering processes which give a different dependence on the coupling constant. 
Generically, two body scattering processes in which a soft bosonic excitation is exchanged 
yield a parametric dependence at leading order as r p ,7 p oc g 2 Tlog(A max /A min ), whereas 
processes in which a soft fermionic excitation is exchanged yield a parametric dependence 
T p , 7 p oc g A T 2 log(q c /gT) /p. The cutoff q c is a scale separating semihard and hard momentum 
transfers, restricted by gT <C q c *C T but otherwise arbitrary, and A max can be chosen of 
order gT. As will be explicitly showed below, the infrared sensitivity to the lower cutoff 
A m in ~ g 2 T entirely disappears from the transport coefficients to be computed. 

On the other hand, the temperature Green's functions for the soft bosonic and fermionic 
excitations are determined by the single particle spectral densities in the hard thermal loop 
approximation, *p L T (u,q) and *A±(u;, q), respectively JR], [IB], |2T|| . These are presented in 



the appendices |B| and |C|. Already, let us note here that only the Landau damping piece of 
these will contribute to the screening of distant collisions. 

B. Kubo formulas and the ladder approximation 

Our starting point is the Kubo formula expressing a given transport coefficient in terms 
of the low frequency, zero momentum limit of the spectral density for the corresponding 
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correlation function. For the electrical conductivity and the shear viscosity the formulas 
are §,0 

1 d 

° = ehS^'^' (10) 

where, as usual, the spectral densities are related to the Fourier transform of the retarded 
correlators by 

/oo r- 
dt / tfW^ x ([J,(t,x), J k (0)])9(t)6 ik , (12) 
-oo J 

/OO P 
dt / d 3 a;e <wt -* q - x ([7ry(t,x),7r^(O)])0(t), (13) 
-oo J 

and the averages are evaluated in the equilibrium grand canonical ensemble. An efficient 
way to compute a retarded correlator IFJ\ (k>, q) (with A denoting collectively the indices 
of the appropriate current) is to exploit the spectral representation for complex frequency 
z. This provides a direct connection with the temperature Green's function HAAfam q)> 
via analytic continuation iv n — > u; + i0 + . Thus, a first step in our basic task is to evaluate 
n^A^n; 0) within the imaginary time formalism. 

After the laborious diagrammatic analysis explicitly performed for the scalar theory J7j, 
the conclusion is that, in order to account for all leading-order contributions to the shear 
viscosity, a set of uncrossed ladder diagrams must be summed. On the other hand, for the 



processes of photon production from a QCD plasma, the authors of Refs. |TT| , [T2| have de- 
veloped detailed power counting arguments which enforce the resummation of the uncrossed 
ladder graphs made of gauge boson rungs. 

A key point for understand the equal footing of this class of diagrams is the presence of 
pairs of propagators carrying nearly the same momenta, which leads to a dependence l/r p 
(or 1/7 P ) by each such a propagator pair. Hence, it is clear that a (n+ l)-loop diagram with 
n uncrossed rungs leads to a dependence proportional to a 2 (l/r p ) n+1 (g 2 ) n , where the first 
a 2 comes from the two external insertions, and each rung introduces a factor g 2 . The derived 
result in both cases [0 and [0] is a linear integral equation for an effective vertex function, 
which is completely equivalent to the linearized transport equation for the problem. Here, 
we will proceed by assuming the dominance of the same set of ladder diagrams, relying on 
a posteriori check of its consistency. 
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A A (iui m + iv n ,iu m ;p) 



w n , 




FIG. 1: Labels of the effective vertex. 



At this point, let us then to introduce an amputated effective vertex Aji(iuj m +iu n , iu m ; p), 
associated with two hard external bosonic or fermionic lines (iu) m + iv n ,p) , (iu m ,p) and 
with an insertion of zero external momentum but non-zero frequency (iv n , 0), corresponding 
to the appropriate current (J, or ny). This effective vertex (Fig. |l|) is the sum of all vertices, 
each one of them with a number n of rungs associated with the exchanged excitations, and 
presumably will encode all collision effects at leading order. The vertex having zero rungs, 
A^(p), does not depend on the frequencies but can depend on the momentum p. 

C. Summation over Matsubara frequencies and analytic continuation 

Let us now to examine the summation over Matsubara frequencies, which generically is 
involved in the evaluation of a temperature correlator H^(iu n , 0), 



where the Matsubara propagators have spectral densities of the form (|) or ([/]). 

To one-loop order, the effective vertex reduces to and the above sum over frequencies 
becomes 



with u n even in any case. Now, the function to be summed is a product of Green's functions 
and one may proceed by expressing the Green's functions in terms of their spectral represen- 
tations given by Eq. ([!]). With this replacement, the resulting expression involves a double 
frequency integral of the product resulting of the spectral densities and the elementary sum 



T^G(iu m + iu n ,p)A A (iu m + iu n ,iu m ]p)G(iu m ,p)Ay(p), 



(14) 




(15) 



1 



1 



(16) 



even, odd m 



7 



where n&j denote boson or fermion occupation factors. However, it is more convenient to 
use an alternative procedure based on contour integration in the z-plane of an appropriate 
function. 

This procedure was used, a long time ago, by Holstein [13] for to derive diagram- 



matically the transport properties of an electron-phonon gas. Here, we closely follow 
the treatment of Holstein. To perform the summation (|ToD, we consider the function 
G(z + ivn, p)G(z, p)ribj(z), and a contour integration Co made of three circuits enclos- 
ing the poles of n b j{z) at the imaginary axis but avoiding the other possible poles of GG 
(in this case there are absent), and the two branch cuts at Imz = and Im(2; + iv n ) = 0. 
This contour may be deformed to C and V as shown in Fig. |2j. The contributions from the 
large arcs T vanish and we are left with the integrals along C. Then, Eq. (|l|) becomes after 
summation 

H(iv n ,p) = ±— I dzG(z + iv n ,p)G(z,p)n b f(z) 

2,7X1 JC 

= T jT H"6,/(0 {[G(iu n + ^p) +G(-iu n + ^p)]G^(^p) 

- G rct (£, p) [G(iu n + e, P) + G{-iv n + £, p)]} (17) 

where £ is a variable specifying the position at the branch cuts. More generally, if the 
function F(z) to be summed has poles in the complex z-plane, the summation formula 
which will be extensively used in what follows is 

ri bj yzi)i\vsyr,z = Zi) zsz ^ / 

even, odd m poles 



-^n 6j/ (0DiscF. (18) 

-co 2,1X1 



Application of this formula to the summation in Eq. (|TJ]) requires the determination of 
the singularities of Aa(z + iv n ,z;p). We argue that the only singularities are two branch 
cuts at the lines Im(z + iu n ) = and lm(z) = 0. This is a consequence of the recurrence 
relation for the vertex with n rungs, 

A^iiUm + iis n , iuj m ; p) = T V / G(iu m + iv n + iu q , p + q) 

v q J ( 27r ) 

xA^ -1) (m; m + iv n + zvg, iu m + iv q \ p + q)G(iuj m + iv q , p + q) ^ ® . , (19) 

g u — %v q 

where p(Q) is the spectral density corresponding to the added rung and (iv q ,q) label the 
momenta running through the loop. For n — 1, it is clear that Eq. (pl|) implies that 
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FIG. 2: Integration contour for the sum T^2 m F(iuj m ). The original contour Co consists of three 
dashed circuits encircling the Matsubara frequencies. The deformed contour C goes along both 
sides of the branch cuts. A possible pole of F{z) at the real axis is depicted for its relevance in the 
summation of the vertex equation. 

has only the singularities of the product of the Green's functions. Then, it follows from 
mathematical induction that A^, and hence A^, inherits the same property. 

Now, we are ready to perform the summation in Eq. fll4|) and the subsequent analytic 
continuation. Making use of Eq. ( |18|) we may write 

-oo Z7TI 

x { [G(£ + iv n )A A {£ + £ - «0 + ) + G(f - iu n )A A (Z - *0 + , £ - iv n )] G adv (0 

- G rct (0 [A A {£ + iu n , £ + tO + )G(£ + iv n ) + A A (£ + *0 + , £ - iu n )G(£ - iu n )] } , (20) 

where the dependence on p is not explicitly exhibited. Next, the analytical continuation 
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%v n — > uj + iO + of the previous expression yields 

T jT || {n 6)/ (£ + ^)G adv (£ + w, p)A^ + w - *0 + , £ - *0 + ; p)G adv (£, p) 

-n 6i/ (OG ret (e + p)A A (£ + w + *0 + , f + *0 + ; p)G ret (£, p) 
- [n 6>/ (£ + w) - n bi/ (0] G rct (e + p)A A (£ + w + *0 + , £ - *0+; p)G adv (£, p)} , (21) 

where we have rearranged some terms by a shift of the integration variable. At this point, 
a simplification arises since the integrand of the above expression contains a large term 
coming from the product G ret G adv . This is due to the fact that the two pairs of poles of this 
product are located at both sides on real axis in the £-plane. Thus, in the limit uj — > 0, the 
contribution of G ret G adv to the integral is inversely proportional to the distance between the 
poles given by the thermal width. The other products G vct G rct and (7 adv (7 adv make a much 
smaller contribution, due to the cancellation between the residues at the poles. 

Therefore, noting that the effective vertex A.a{£, + i0 + , £ — i0 + ; p) is a real quantity 1 , and 
n 'bj(0 = ~P n bji.0(^-^ zn bj(0)^ the required zero frequency slope of the retarded correlator 
may be written as 

^ 3 P A (o)/ N f°° 
xG rct (£, P )A A (£ + *0+, £ - z0 + ; p)G adv (£, p), (22) 



£imn^,o; 



^ = ^/g S AS)(p)£| nw (0(l ± n w (0) 

adv / 



where the prefactor ( accounts for the symmetric factor associated with the one-loop dia- 
gram corresponding to IL^. This factor is 1/2 or 1 depending on whether the G functions 
correspond to a self-conjugate field or not. It is important to notice the correspondence 
between Eq. (|22"D and the expression for a transport coefficient in terms of 5riA,±{p), de- 
noting the departure from equilibrium of the single particle (anti-particle) density function. 
Anticipating a contribution of G Tet G adv oc J2± <K£ T |p|)/T p , such a correspondence is clear 
with the appropriate identification 

6n A ,±(p) oc ^A A (± |p| + *0 + , ± |p| - i0+; p). (23) 
1 P 



1 This follows from mathematical induction, taking into account that the zero-order vertex AjP (p) is real. 
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III. THE TRANSPORT EQUATION FOR SHEAR VISCOSITY IN A0 4 THEORY 



With the purpose of making clear the basic procedure to be used in gauge theories, we 
present a simple alternative derivation of the transport equation for the scalar theory in the 
form previously obtained by Jeon [0. The starting point is the linear integral equation for 
the effective vertex before analytic continuation 

o 1 f d^Q 

Aij(iuj m + iu n , iuj m ; p) = A^ 0) (p) + ~ T J2 TT\i G ( iuJ rn + iv n + iv v P + q) 

xAij(iuj m + %v n + iis qj iuj m + iv q \ p + q)G(iuj m + iv q ,p + q) ^} , (24) 

where the factor 1/2 is due to combinatorics, and A^(p) = 2p 2 (p l p ) ' — 8^/3) 2 . This equation 
is shown in Fig. ||. The spectral density p(Q) associated with the bubble in the rung is given 
by 

f . r J 3 h i 1 1 

(25) 



p(g ,q) = 2Im {a 2 T£ j 



ii/ q =q°+iO+ 



(2n)3(v g + v k )i + \k + < 1 \ 2 vl + V 

By performing the previous sum (see appendix |A|), and using the expression for the free 
Wightman function, 

G>(P) = (l + n b (p )) 2vrsgn(p°) 6(p% - P 2 ), (26) 
one may write the spectral density for the rung in a symmetric form 

" ra = ^/H G » >( - - A ' )G » (A '»- (27) 

The relation of this spectral density with the kernel ^Boitz in the treatment of Jeon [^] is 
therefore 

P(Q) = ^yWz(-g). (28) 

Note that p(Q) is odd in q°. 

The formula (|Tj^) can be applied in order to perform the summation in the vertex equa- 
tion (p4j). The term associated with the discontinuities gives (omitting the prefactor 1/2 



2 The factor 2 is due to the sum over the two permutations of the momenta corresponding of the scalar 
field. 
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FIG. 3: Equation for the effective vertex in the A</> 4 theory. The double lines represent the propa- 
gators of the hard particles including the 0(A 2 ) thermal width. 



and p{Q)) 



30 <c\ I ^n,P + q) 

77- n b {0 < 5 : — 

00 2m I 4 — g u — iuu m — iv n 



x 



G adv (£, p + q)A^(£ - *0 + , i - iv n ) - G rct (£, p + q)A^(£ + *0 + , £ - 
C(C + Wn,P + q) 



+ 



X 



G adv (£, P + q)Atf(£ + ii/ n , £ - *0 + ) - G rct (£, P + q)A^(e + ^n, £ + *0 + )] } , (29) 

and the pole term associated with the spectral representation of the rung gives 

n b (q°)G(q° + iu m + iu n , p + q)Ay (g° + iu m + iv n , q° + iuJ m )G(q° + iu m , p + q), (30) 

where we have omitted the dependence on p + q in the vertex. After this summation, the 
analytic continuation iuj m + iv n — > p° + u + i0 + , iu m — > p° — i0 + can be explicitly performed 

A^-iO+e-M'n) -> A^-*0 + ,£-cu-*0 + ), 

Ay(£ + i0+£-ii/ n ) -> A^(£ + i0 + ,£-u;-;0 + ), 

+ -i0 + ) -> A Ji (e + ^ + «0 + ,e-«0 + ), 

Ay(f + iz/„,£ + i0 + ) -> A ii (cj + £ + i0 + ,£ + i0 + ), 

G^-u^p + q) - G adv (£-cu,p + q), 

C(e + ^„,p + q) - G rct (£ + cu,p + q). 

If we neglect the products G adv G adv and G rct G' rct as argued before, the discontinuity contri- 
bution becomes 



J-00 

MO ^fc(£ + ^) 

£ - p° - g° + i0+ £ - p° - q° - i0+ 



(31) 
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and the required limit u — > 0, after the £ integration, reduces to 

- n b {p° + q°)G TCt {p° + q°,p + q)A^(p° + q° + t0 + ,p° + q° - iO + )G adv {p° + q°,p + q). (32) 

Adding the pole contribution, one finds the integral equation satisfied by the effective vertex 

Va{P) = AS?(P) + 11 ^ G-\P + Q)V tJ (P + Q)G^(P + Q) 

xp(Q){n b (q°)-n b (p + q )], (33) 

where we have defined TV,(p°, p) = Aij(p° + i0 + ,p° — i0 + ; p). 

To present a more explicit form of the transport equation, it remains to analyze the 
product G ret (P + Q)G adv (P + Q). If the momentum P (or P + Q) is nearly on-shell, the 
contribution associated with the sum of the crossed products reduces to a product of two 
Lorentzians peaked at different values or p°. In the limit A — > 0, this fact enables us to 
neglecting the piece G^G^ + G^G^, so we may write 

G rct (P + Q)G adv (P + Q) = V - 1 , 2 — n ^-—^ — ,, 2 ^ 2 (34) 

± 4| P + q| 2 (p° + q° t |p + q|) 2 + r 2 +q 

which, for A — > 0, behaves as 

G ret (P + Q)G adv (P + Q) = 71 (5(p° + q° - |p + q|) + 5(p° + q° + |p + q|)) 

4|p + q| 1 p+ q 

-^rk^ 2l ^ (p ° +9 ° )2 - |p+q|2 )' (35) 

Since the imaginary part of the retarded self-energy is odd in the frequency, the thermal 
width can be replaced by — sgn(p°+g°) Im E ret (P+<5)/(2 |p + q|), and Eq. (RH) then becomes 



G^(P + Q)G^(P + Q) = - 2ImSrct 1 (p - 27rsgn(p° + q°)6 ((p° + g ) 2 - |p + q| 2 

= ~ (1 + n h (jP + q°)) 2 ImS"*(P + Q) G ° >(P + Q) ' (36) 

With the aid of the identity 

n b (q°) - n 6 (p° + q°) = n b (q°) (l + n b (p° + q )) (l - e"^) , (37) 
one arrives to the final form of the transport equation 

v ii{ P) = 4 0) (p) - (i - *-**) J H W-Q) + Q) 2 Z'^ P Q + Q) ' < 38 > 
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in complete agreement 3 with the results of Ref. 0, where the equivalence to the linearized 
Boltzmann equation has been proven. The function T>ij(P) is real, even in p°, provided that 
p(Q) is odd in q°. 

The last step in the evaluation of shear viscosity is the computation of the integral 
for which only the on-shell effective vertex T>ij(p° = p, p) is required. The insertion of (3"B) 



with P + Q replaced by P into and the multiplication by 1/10 yields the shear viscosity 
compactly written as 

IV. SOFT CONTRIBUTIONS TO THE TRANSPORT EQUATIONS IN GAUGE 
THEORIES 

In the derivation of the transport equations that we try, we want emphasize the role 
played by the screening effects in the regularization of some infrared divergences which arise 
in the small momentum transfer region. The importance of the dynamical screening in 
transport phenomena in gauge theories was first recognized in Ref. |]], where a linearized 
collision integral free of infrared divergences was stated by using only screened interactions 
mediated by gauge bosons. Our aim here is to derive transport equations by means of the 
summation of a restricted set of ladder diagrams, which includes only a specific type of 
rungs. These rungs will consist of appropriate effective bosonic or fermionic propagators in 
the hard thermal loop approximation, thus accounting for the screening of distant collisions 
between all different plasma constituents. Obviously, this approximation does not include 
the effects due to close collisions and its use requires the imposition of an upper cutoff q c 
in the integrals over the momentum transfer ||15||. As a check of consistency, we will verify 



that the coefficients of the UV-logarithmic sensitivities to q c match to the IR-logarithmic 
divergences previously computed by Arnold, Moore and Yaffe {§] in their treatment of the 
hard contributions to the collision integrals of the Boltzmann equation. 

To elucidate more precisely what type of rungs dominate in the ladder summation which 
gives the logarithmic accuracy of transport coefficients, it is necessary to estimate the power 



Our notation docs not agree with Q|, but the conversion is direct: AV°' — * 21^,1)^ — » 22?^ and Im57 e 



L (0) 
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counting size of the spectral densities which, in our treatment, are associated with the ladder 
rungs. For instance, consider the spectral density of the rung in the scalar theory of the 
previous section, Eq. fl2~T|). Clearly, it is 0(A 2 ) from the two vertices, and the contribution 
0(A -2 ) from the product G ret G adv in Eq. ( |5BD causes the integral term in the vertex equation 
to have a net g° contribution, as the zero order vertex. 

In a gauge theory, there are two possible topologies with the same power counting which 
potentially contribute to the leading logarithmic order 4 . They are illustrated in Fig. |]. The 
spectral density of the soft gauge boson exchange ladder, labeled as (a), is 0(1) because 
of the g 2 suppression from the two vertices, and a g~ 2 T~ 2 contribution from the spectral 
density of the soft gauge boson. In the Q integration of the vertex equation, similar to 
Eq. fl33|), there is a factor n^q ) = 0(g~ l ) and a factor 0(g~ 2 ) from the imaginary part of 
the self-energy in the denominator of the product A ret A adv . Finally there is a g 3 suppression 
from the soft integration c? 3 q (the contribution of the integration dq° is cancelled by a delta 
term from A ret A adv , see Eqs. (^) and (|5TD below). Hence, the net contribution of the 
integral term in the vertex equation is 0(1), the same as the zero order vertex, when the 
spectral density of the rung is also 0(1). 

Similarly, when the horizontal internal lines are soft gauge boson propagators, the spectral 
density of the box type ladder rung, labeled as (b) in Fig. ^, is 0(1). To perform this 
estimation, it is helpful to write the general form of the spectral densities in terms of the 
squares of the 2 <-> 2 scattering amplitudes |.M| 2 . Such formulas are similar to Eq. ( |2"T| ) 

Pa (P',P) oc . * Q . f d 4 Kd 4 K'5^(P + K - P' - K')\M\ 2 A > (-K')A>(K), (40) 
Pb (K',P) oc 1 / d 4 Kd 4 P'5^(P + K-P'- K')\M\ 2 A>(-P')A>(K), (41) 

where P, K, P', K' are the on-shell hard momenta for the particle entering into the scattering 
processes, and the labels for momenta have been chosen with the aim to use the same 
\A4(P, K, P', K')\ 2 within the integrations 5 . For soft momentum transfer Q = P' — P = 
K — K', this reduces to 

\M(P,K,P',K>)\ 2 oc gVk 2 \*A L (Q) + (p x q) • (k x q)*A T (Q)| 2 ~ O(g ) , (42) 

4 I thank the referee by pointing this out. 

5 We closely follow the notation of Ref . || . 
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(a) 



K I 



t P' 



(b) 



FIG. 4: Two ladder rung topologies of the same power counting. The dashed lines correspond to 
the cuts giving the spectral densities p a (P',P) and pb(K',P). 

where *A Lj t denotes the longitudinal or transverse retarded boson gauge propagator. Note 
that P' and K are loop momenta associated with the integrations in the vertex equation, 
and P corresponds to the fixed momentum of the effective vertex. In this language, it is 
easy to recover the power counting size of p a (P', P): 

Pa (P',P) oc q° [°°dkri bJ (k) /^ k 5(g°-k.q)|^| 2 oc^VT 2 ^[|A L (g)| 2 



1 



2\ 2 
% 



i^t(q)i 



2 2 

cx g p 



% 2 ' 



0l(q) + i - ^ Pi{Q) 



0(/).(43) 



For the estimation of p\>(K',P), the denominator contributes as g° even though k'° and p° 
are hard energies. The reason for this is that, in general, their difference does not correspond 
to a soft exchange of energy, as shown in Fig. f|(b). Consequently, ph(K' , P) = 0(1) also. 

However, in the case of the shear viscosity, the subsequent integration in the vertex equa- 
tion over the directions of k' and the q° integration cancels the contribution of this box type 
contribution. To understand this cancellation, we note that for small momentum transfer 
q p,k, the energy 6 function which remains after the three-momentum p' integration in 
Eq. ( [4"ID may be expressed as 5(q° — k' • q), and 



M(P,K,P',K') oc g 2 pk 



*A L (Q) + *A T (Q) 1-^ cos 



(44) 



where <fi is the angle between the planes (p, q) and (k', q). On the other hand, the angular 
dependence of A 7r (K / ) within the integration over K' in the vertex equation is given by the 
contraction of A n (K') with A n (P), which is proportional to P2(cos8 p k>), with P2 the second 
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Legendre polynomial. For small q, the required cosine becomes 

cosV-^ + (^--f)™^- ( 45 ) 

Without screening effects, the boson propagators reduce to *Al(Q) = — l/q 2 and *At(Q) = 
— l/(q 2 — q 2 )- Thus, the angular integration over the azimuthal angle of k' followed by the 
k° (or q°) integration turns out to be proportional to 

P dq° r d<pP 2 {cos6 pk ,){l-cos<p) 2 = 0. (46) 

J-q JO 

Here, the integration over the spatial range — q < q° < q is enforced by the 5(q° — k' • q). 
This means that the unscreened box rungs do not contribute to the leading logarithmic 
order and, consequently, box ladder rungs made of two soft gauge boson propagators are 
irrelevant in order to compute the soft contribution to the shear viscosity. 

For the electrical conductivity at zero chemical potentials, the cancellation of these box 
ladder rungs comes from the charge-conjugation invariance. For each fixed P, we have two 
contributions to the integral term of the vertex equation, corresponding to the insertions of 
fermions and antifermions through the loop. These contributions are of opposite sign, so 
their sum vanishes. 

Similarly, there are other box ladder rungs in which the horizontal lines are soft fermion 
propagators, which are comparable to the soft fermion exchange ladders. They do not 
contribute either to the leading logarithmic order because of the cancellation under angular 
integration of their unscreened counterparts. Now, the unscreened squared amplitudes are 
| M. | 2 oc (1 — cos <p)/q 2 Q , and the relevant integration is 

f 1 dq° f 277 d(j)Pi (cos 6 pk ,){l -cos <p) =0, (47) 

J-q JO 

with 1 — 1,2 for the electrical conductivity and the shear viscosity, respectively. 

Next, let us consider the equations for the effective vertex of a given fermionic species 
and the effective vertex h 3 ™ for a given boson gauge. Here, the superindices (jm) denote 
spatial indices corresponding to the boson propagators in Coulomb gauge to be joined to 
the vertex, and A collectively denotes the indices corresponding to the insertion of ji or TTjfc. 
We still does not explicitly indicate any spin, color or flavor indices. The equations for the 
effective vertices are illustrated in Figs. | and These equations sum all ladders consisting 
of rungs made of a HTL gauge boson propagator, and also a HTL fermion propagator. The 
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p\ 



p\ 



X 



FIG. 5: Equation for the soft contribution to the fermionic effective vertex. For the case of the 
electrical conductivity to leading order in e, the soft boson exchanged corresponds to a photon and 
the last graph is zero. 



FIG. 6: Equation for the soft contribution to the bosonic effective vertex. In the case of the 
electrical conductivity, all graphs are null. 

exchange of a soft gauge boson takes place in scattering between fermions and gauge bosons 
in non-abelian theories, while the exchange of a soft fermion enters in an annihilation process 
into a fermionic pair, and the inverse process of creation. 

By following a completely similar treatment to that we have presented for the scalar 
theory for the summation and analytical continuation to real frequencies, one arrives to the 




+ 
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coupled equations for the effective vertices 



A A (P) = A^(p) + g 2 



I 



d 4 Q 



^S ret (P + Q)A A (P + Q)S adv (P + Q)y 



(27T) 4 



+ 9 




x*p(c,q) 7 " nf ( q °) + n b (p° + q ) , 



(48) 



and 



A7(P) = A 




Q)V™ n {-P,-Q,P + Q) 




Here *p^, v (Q) and *p(Q) denote the spectral densities of the soft gauge boson and the fermion 
and, as usual, A^(P) = Aa(p° + iO + ,p° — iO + ;p). Note that we have used the parity 
properties *p(Q) = *p(—Q) and q°) — nf(—p° — q°) = — [n/(g°) — n/(p° + q )]. These 
equations are entirely similar to the transport equation for the scalar theory. As there, the 
occupation numbers fibj(p° + q°) arise from the branch cut contributions, and the identities 
n b,f(£ + niti) = — Ji/,&(£) with n odd have been used, if necessary. Further progress requires 
to examine in detail the explicit structure of the products 5 ,ret A^5' adv and G ret A J 4G adv . 

The substitution of the two pieces of the Green's functions in the product S TCt S adv or 
G ret G adv yields four terms, two of them A^*A^ V or G^G^ can be directly dropped in 
the limit g — > 0. Now, by computing the soft contribution to the transport properties, we 

have only to consider one of the two terms, ++ or , depending on whether the external 

momentum corresponds to p° = p or p° = —p. This important simplification is due to the 
fact that, when the external momentum P is on-shell, let's say p° = p, the other sheet of the 
mass-shell can not be connected by a soft momentum transfer. Clearly, a hard momentum Q 
is required in order to create two propagating on-shell particles both with hard momentum. 
In contrast, for the case of the scalar theory, we have retained both terms in Eq. ([35|) because 
there, we have not made distinction between soft and hard momentum transfers. Thus, when 
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P is on-shell and Q <C P, we can make the approximations 



Af(P + Q)A* dv (P + Q) = — 5(p° + <?° =F |p + q|) — " 



7p+q 



7p+q 



%°=Fp-q), (50) 



and 



Gf(P + Q)Gf v (P + Q) ~ 



7T 



%°Tp-q). 



(51) 



4p|p + q|r p+q 

On the other hand, in the high temperature limit when all masses are negligible, the the- 
ory is chirally invariant. As a consequence, the projection operators h±(p) can be expressed 
in terms of simultaneous eigenspinors of chirality and helicity, 



Mp) = u(p,+)u(p,+) + u(p,-)u(p,-), 

Mp) = «(p>+Mp>+)+«(p,-)«(p-)> 



(52) 
(53) 



where the u (v )-spinors have the chirality equal (opposite) to the helicity. For example, in 
the chiral representation, with the momentum along the 3-axis 



w(p,+) 




1 

V / 



, «(p,-) 



1 
o 

V / 






V / 



, u(p,-) 






V 1 / 



(54) 



In the case of the electrical conductivity, the zero-order effective vertices are 



A 



(0)jm 



0, 



(55) 
(56) 



where q s is the charge of the fermionic constituent s in units of e. For the shear viscosity, 
the insertion at zero momentum of a spatially transverse energy-momentum tensor yields 
the zero-order vertices 



A S?(P) = ^(tV +7V-^7-P^), 



A 



(0)jm 
ik 1 



p) = 2[p l p k -^-5 ik )5 3 



(57) 
(58) 



Both of the fermionic vertices are linear in the 7-matrices, and this linear dependence is 
preserved by summing ladder diagrams because each added rung does not introduce any 
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extra dependence. Thus, the fermionic effective vertices, which appear sandwiched be- 
tween the projection operators h±(p) and h±(p + q) require to consider a 7-matrix be- 
tween eigenspinors of all possible chiralities and helicities. Obviously, the combinations 
between spinors of different chirality, as u(±)7 a 'm(=p) or v(±) / j fi u(±) , are zero. Moreover, 
the combinations mixing the particle and antiparticle mass shells, as v(p + q, ±)7 M w(p, =f) 
or u(jp + q, ±)7 /i v(p, =)=), do not need to be retained since they can not be connected by a 
soft momentum transfer. Therefore, we are left with the chirality-independent combinations 

U(pTq, ±hVF>, ±) = 1 + 0(q/p), «( P Tq, ±hMp, ±) = f + 0(q/p), (59) 
v(pTq, ±)l°v(p, ±) = 1 + 0(q/p), tJ( P Tq, ±)tMp, ±) = ~f + 0{q/p), (60) 

corresponding to the lepton (quark) and the anti-lepton (anti-quark), respectively. 



A. Soft contribution to the transport equation for the electrical conductivity 

Now, we are ready to write more explicitly the transport equations ( f48|) and (fi9|). For 
the case of the electrical conductivity, it is suggested to define the non-amputated on-shell 
vertices for a given charged fermionic species s 

Vf(p) = u( P) ±Mp°=p,pMp,±), (61) 
vf{p) ee U(p, ±)At(p° = -p, -p) v (p, ±), (62) 

with the corresponding zero-order vertices which follow from Eqs. (|59D and (|60D. Here, the 
factor g 2 in the front of the integral in Eq. (|48|) must be replaced by q 2 e 2 . At zero electrical 
charge, the bosonic effective vertex for electrical conductivity does not enter because its 
zero-order vanishes, and Furry's theorem ensures the vanishing of the term within the trace 
in Eq. (|49|). This leaves a single decoupled equation for the fermion vertex. 

Next, we fix the external frequency p° = p. After multiply both sides of the equation ( ^8|) 
by the u-eigenspinors, and expand the integrand using (|52|)and (|50"D with ++, we obtain a 
term proportional to 

8{q° - p ■ q)u(p, ±)7"«(pTq, ±) u(pTq, ±)7^(p, ±) * P^(Q) = 



S(q° 



p - q) 



0[q/p], (63) 
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where the delta function has enforced a space-like momentum transfer momentum corre- 
sponding to Landau damping, and we have used the approximations (|59|). Finally, the 
substitution A l (|p| + p • q, p + q) = A*(|p + q| , p + q), valid at soft momentum gives the 



equation 



£>f (p) 



q s ep l + q 2 s e 2 



d 4 Q 
(2tt)* 



7r5(q° 



p-q) (-0 - \ + n f{v) 



X 



fa(Q) + 1 - 2a (3 T (Q) 



te + q) 



7 P +q 



(64) 



where we have neglected the 0(q°/p) terms in the occupation numbers. Using the parity 
properties of the zero-order vertex and noting that the spectral density *p M ^(Q) is odd in q°, 
one may easily show from (f49|) that T>\ (p) = — P| + (p). This equation is formally similar 
to the integral equation for the complete leading order of the photon emission rate from the 



quark-gluon plasma, which has been derived in Ref. JTT 



Following a similar treatment to that used by Arnold, Moore and Yaffe [FT], we may 
reduce a bit more the transport equation. This proceeds by insertion of the integral for the 
thermal width of a hard fermion. As we have already said, the thermal width receives the 
leading order contribution from two-body scattering between fermions by exchange of a soft 
photon. For a given fermionic species, it is given by the integral ||T6| 



T ( 2 ), 
'P 



2 2 



d*Q 
(2vr) 4 



ir8(q° 



p ■ qj 



T 



Pl(Q) 



i - 



q 2 



Pt(Q) 



(65) 



Note that the value of this integral remains unaltered when the term T/q° of the integrand 
is replaced by T/q° — 1/2 + n/(p). This is due to the fact that the added term — l/2 + n/(p) 
is here multiplied by an odd function of q° and, consequently, is irrelevant. However, this is 
not true for the corresponding term in (0), since the vectorial dependence of Pj on p + q 
can give odd contributions in p • q. 

Besides the dominant contribution, the thermal width receives another contribution com- 
ing from two-body conversion processes which also give rise to a leading-log term in the 
transport coefficients. This contribution to the thermal width corresponds to the imaginary 
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(b) 



FIG. 7: Subleading contributions to the self-energy of a hard quark (a) and a gluon (b). 



part of the diagram of Fig. ^, and it is expressed by the integral 



7- 



(4), s 



X 



2pJ (2v 

/3 + (Q) (l-^)+/3_(Q) U + £ 



8np 



;i+2n 6 (|p|)) 



In 



\/2< 



- 1 + ln2 



(66) 



where g c is the upper limit of integration which separates semihard and hard momentum 
transfers, and u 2 s = q 2 e 2 T 2 /8 is the plasma frequency for the fermionic species s. 
By inserting the expression ( j65|) into (|64D , after defining the quantities 



Xl (P) 

one obtains a more convenient expression 



qyp' 



^r + ( P ) 

7n 



TTd{q -p-q) 



(67) 



(2vr 



% 2 



/3 L (Q) + 1 - 2L\ (3 T (Q) 



2 

xf (p) -xf (p + q) 



n/(p) 



(68) 



for the soft contribution to the transport equation in the case of electrical conductivity. 

An important point to notice here concerns to the absence of some subleading corrections 
to the hard particle width, whose size is comparable to that included in Eq. (IBB]). Such 
corrections of 0(g 3 ) and 0(g 4 In g~ l ) arise from the subleading part of the bosonic occupation 
number n&(g°) in the integrand of Eq. (|5"5D , and they are associated with 

(69) 



T 1 00 q u 



To understand this absence, consider the transport equation ( |6"31) when T>/ r y has been re- 
placed by x, as in Eq. (0), and the kernel still contains n b (q°) + ri/(p). Now, if in the 
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product 7 P x(p) one were to include the subleading terms from Eq. (0), the same terms 
would have to be included into the occupation number in the kernel of the integral. For 
small q, there is a cancellation between them, and the resulting transport equation would be 
the same as that in Eq. with T/q° — 1/2 + rif(p) replaced by nb(q°) + iif(p). However, 
when the momentum p is hard, n/(p) = 0(1), and other pole terms different from zero in 
the occupation number are 0(g). Hence, they may be ignored. 

B. Soft contribution to the transport equations for shear viscosity 

Next, we consider the shear viscosity in a SU(N C ) gauge theory with Nf fermion fields in 
a given irreducible representation r of dimension d(r). The generator matrices are denoted 
by t%, and the normalization for this representation is defined by the constant C(r) in 
tr(t^t b r ) = C(r)5 ab . The quadratic Casimir operator is denoted by C^r). 

Since the gluon effective vertex A^ n,ab is always joined to a pair of transverse projectors 
of the gauge propagators, it is useful to define the non-amputated on-shell gluon vertex T>\ k 
by 

p^'(p)Ar a V = IpI , p)Pt s (p) = VUp)Pi s (p)S ab , (70) 

and the zero-order vertex corresponding to Eq. (|58|) , Z>^ 9 (p) = 2(p l p h — 5 lk p 2 /3). On the 
other hand, and as before, we define the the non-amputated on-shell quark vertex T>\ k by 

V\ h (p) = u(p,±)A ik (p° = p,p)u(p,±), (71) 

and the zero-order vertex corresponding to Eq. (|57D, (p) = \p\(p l p k — 5 lh /3). We 
note that the similar vertex for anti-quarks, Tfi = vA(—P)v, turns out to be the same as 
for quarks, as it is easily checked by examining the corresponding zero-order vertices, and 
noting that each added rung in the iteration does not break this property. 

Now, in order to derive the transport equations for T>f k 9 ', it only remains to perform the 
appropriate contraction of Eqs. ( [ISP and ( f49"D with a pair of w-spinors and a pair of transverse 
projectors, respectively. Using the approximations 

u(p, ±) 7 m 0+(Q)Mq) + P-{Q)h-{q)] IMP, ±)Pt U (pTq) = 

/3+{Q) (1 - P • q) + /3_(Q) (1 + p • q) + 0(q/p), (72) 

tr { 7 m MpTq) 7 n [P+(Q)h+(q) + [3-(Q)h_( q )]} P™(p)P™(p) = 

\j3+(Q) (1 - p • q) + (3-{Q) (1 + P ■ q)] P£(P) + 0(q/p), (73) 
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and 



f^f^v^iP, g, -p - Q)v mm (-p, -g, p + q) v^(g)p^( P Tq)P^(p)pr(p) = 

ViV c p 2 [ft,(g) + (l - p ■ q 2 ) /? T (Q)] P-(p)5 ab + 0(g/p), (74) 

one finds 

VUp) = V^(p) + g*C 2 (r) [ *Q n5( q ° - p • q) B(Q) f £ - ~ + n f (p)) ^MR±A 

J (2vr) 4 \q u 2 y 7 p+q 



and a completely similar for the gluon vertex 



2%(P) = 2>£ } '(p) + 9 2 K | P | / vr5(g° - p • q) B(g) f-^--- n b (p) 



(2^-- - — , Vg0 2 j | p + q |r p+q 

q) 

(2^"^ - — — V2 '^7 7p+q 



+2, 2 iV / C(r) / ^(g - p - q) ^(g) (I - n/ (p)) ^ P + q) , (76) 



where 



s(g) = /? L (g) + (i " ^) /W) , (77) 

.F(g) = /? + (g) (i - 7) + /MG) ( x + 7) • ( 78 ) 

Note that we have approximated the occupation numbers by their expansions up to 0(q°) 
corrections. The prefactor 2 in front of the piece g 2 NfC(r) in the last term of (|76| ) is due 
to the two possible orientations for the momentum in the quark loop. The remainder of the 
previous equations is almost obvious, after performing the approximations ( |50l) and (|5l|) . 

Finally, as for the case of electrical conductivity, in order to make more clear the closed 
resemblance of these equations with Boltzmann-type equations, we may define the quantities 
which will turn to correspond to the deviations from the equilibrium distribution functions 
of quarks and gluons, 

X? fc (p) = ^, (79) 
7 P 

X«(P) = S#- (8°) 
Z IPI 1 p 

Then, the substitution of thermal widths by their integral representations yields the final 
form of the soft contribution to the coupled transport equations for these quantities. They 
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become 



\v\m k -\* ik ) 



iC 2 [r) J n5(q° - p • q) B{Q) (^-\ + n/(p)) bd(p) - X? fc (p + q)] 



(2tt 

S 2 C 2 (r) /l , , A f d 4 Q 



2|p| V2 + «6(P)) / ^ ^(9° " P ■ IxSb(p) " X?,(P + q)] , (81) 



and 



^ c / „6(q° - p • q) f ^ - 1 - n 6 (p)J [ X f fe (p) - xUp + q)] 



g 2 N f C 2 (r) (I , \ f d A Q 



n f ( P )) J -JUttf ~ P ■ <l) IX&(P) " X?*(p + <!)] • (82) 



IpI V2 ^vj (2 

V. EXTRACTING THE LEADING-LOG TRANSPORT COEFFICIENTS FROM 
THE SOFT CONTRIBUTION 

Rotational invariance fixes the form of the x-functions. In the case of electrical conduc- 
tivity, they have the form 

xf(p)=P i X s ~(p), (83) 

and for shear viscosity, 

X%\V) = {p l p k - \t lk )x q ' 9 (p)- (84) 

When |q| -C |p|, we may to approximate the integral collision terms in the transport 
equations by a derivative expansion of the x-functions, hi a similar way to the procedure 
which leads to the derivation of Fokker-Planck equations in the context of classical plas- 
mas |T7|, With the aim to extract the leading- log terms for the transport coefficients, we 
only require the insertion into the transport equations of the second-order approximations 
for the quantities ft 1 xt (p) ~ X\ (p + o) an d p l p k [xtk 9 (p) ~ xik^P + c 0}- At the same ac- 
curacy, terms of the type \xlk\p) - xfk^P + o)] are replaced by [xlkiP) ~ xlkip)]- Triese 
substitutions produce a combination of derivatives of the x-functions, whose coefficients are 
integral expressions involving the Landau damping piece of soft spectral densities. Here, 
the leading-log terms will arise from to the logarithmic dependence on the upper limit q c . 
Although the computation of the complete leading order in g of these integrals can need 
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numerical quadrature, the leading-log order is easily obtained analytically by means of the 
formulas based on sum rules given in the appendix [B| 

After the angular integration of the Eqs. flST]) and (p2|), and the subsequent q° integration 
with the aid of Eqs. ( |B11| ) and ( |B12| ), the logarithmic terms finally combine to give the 
following expressions in the case of shear viscosity 



P 



g 2 C 2 (r)m 2 D T 



16n 

9 2 C 2 (r)u; 2 



Ins" 1 \x q (p)" + 



2 1, 

--- (l-JJn/Cp 
p 1 



x q (p)' - - 2 x q {p) 



inp 

,2 AT ™2 



\ng- 1 (l + 2n b (p))[ X q (p)-X 9 



(85) 
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g N c m D T _J 
— \ng <x g {P) + 
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2 1, 

--- l + 2n 6 p 
p 1 



x 9 (pY - %x 9 {p) 



g2N fJP )UJ2 ° lnrl (1 " 2nf{p)) [ ~ x9{p) + x9{p)] ' 



(86) 



where lng 1 is the term coming from either In (jj^j or In with logarithmic accuracy, 
and 



g 2 T 2 

3 

9 2 T 2 
8 



m 



D 



[N c + N f C(r)\, 
C 2 (r). 



(87) 
(88) 



Quite remarkably, the terms (l±2nb ) /(p))/T multiplying the first derivatives and gener- 
ated by the zero-order in q° of n b (q°) =F n b j(p + <7°)j allow us to write a single functional of 
X 9 {p) and x 9 {p) whose variation leads to above equations for shear viscosity in the leading- 
log approximation. This functional dpC(x, x') turns out to be exactly the same as that 
previously found by Arnold, Moore and Yaffe in their derivation of the leading-log terms 
of the linearized collision integral of the Boltzmann equation. After multiplication of both 
sides of Eqs. (|85|) and (j36|) by 2d(r)Nfp 2 rif(p)(l — rif(p)) and d{G)p 2 n b {p){l + n b {p)) respec- 
tively, one finds 

C( X , x') = -2d{r)N f p 3 n f (p){l - n f {p)) X q {p) - d{G)p z n h (p){\ + n b {p)) X 9 {p) 

g 2 N f uj 2 d(r)C 2 (r) _ x 2 
+ q ln# pn f (p)(l +n b (p)) [x q {p) ~X 9 (P)\ 

+ g 2 N f m 2 D Td(r)C 2 (r) ^ _ ^ r^)* + g ^ )5 



16n 

+ g 2 N c mlTd(G) ln ^-i^ (p)(1 + Mp)) r^)* + 



32vr 

where we have used d(r)C 2 (r) = d(G)C(r) with d(G) = N 2 — 1. 



(89) 
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Using the expressions fl22|) and (|ll|), and noting that (p l p k — 5 tk /3)(p l p k — 5 lk /3) = 2/3, 
we see that the contribution to the shear viscosity of each Dirac fermion is 

3 r°° 

1 dpp 3 n f (p)(l-n f (p)) X q (p), (90) 



15tT 

while each gauge boson contributes as 



8 



oo 

3 



dpp n b (p)(l + n b (p))x 9 (p). (91) 



30vr 2 Jo 

Hence, the shear viscosity is exactly the value of the functional — 3/(15ir 2 ) / °° dp £ for the 
values of x solving the motion equations. 

For the electrical conductivity, the above procedure applied to Eq. (^) yields the expres- 
sion 



2 1, 

---(1-271/0?)) 



+ -^-\ne- 1 (l + 2n b (p)) X s '(p), (92) 

where the Debye mass for the photon and the plasma frequency for the charged fermionic 
species s are given by 

e 2 T 2 



m 



D ~ 3 



Uptons + d{r) ]T q 2 \ (93) 

quarks J 

o 2 e 2 T 2 

^ = (94) 

This equation may be retrieved by varying the functional / °° dpC s (x, x')> where C s may be 
chosen as 

J 2 

C s = -p 2 n f (p)(l -n f (p))x s (p) + j^\ne' 1 pn f (p)(l + n b (p))x s (p) 2 

+ ^lne- 1 n f (p)(l-n f (p))[p 2 x s '(p)' 2 + 2x s '(p) 2 }. (95) 

OZ7T J 

Now, the contribution of each charged Dirac fermion to the electrical conductivity is 

8 f°° 

— dpp 2 n f (p)(l -n f (p))x s (p), (96) 

07T Z JO 

which implies that this transport coefficient is given by — 28/(3tt 2 ) J2s Jo° dp £ s at the sta- 
tionary values of x s ■ 

This completes the derivation of the leading-log terms for the transport coefficients. The 
agreement with the results of Ref. || is complete when the conversions x 9,9 ~2PXamyi 
X s — * —^PXamy are performed, and the sum over charged species is restricted to leptons. 
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VI. CONCLUSION AND PROSPECTS 



In this paper we have shown how to derive transport equations in some relativistic many 
body theories from the summation of uncrossed ladder diagrams within the imaginary time 
formalism. The procedure is similar to the one used in a quite remote past by Holstein [13j 
and, for the case of the scalar field theory, yields the correct results previously derived by 
Jeon J7|. In this treatment, one first identifies the analytic continuation to real frequencies 
for the effective vertex, and then writes the integral equation which it satisfies. The two 
relevant quantities in the vertex equation turn out to be the imaginary part of the self-energy 
and the spectral weight of the ladder rung. 

For the case of gauge theories, we have derived a couple of transport equations (|68|), fl8lD 
and ( |32D by resummation of the ladder series, whose graphs are made of the rungs associated 
with resummed propagators in the HTL approximation. The kernels of the integral equa- 
tions derived in this way are the same as the integrands of thermal widths for hard particles 
to leading and next to leading order. By extracting the logarithmic terms, the transport 
equations turn out to be differential equations similar to Fokker-Planck equations which 
appear as approximations to the collision integrals in Coulomb plasmas. The 0(g 4 In 1/g) 
thermal widths in these equations are seen as damping terms, similar to those which arise 
when the Boltzmann equation is treated in the relaxation time approximation. These dif- 
ferential equations are the same as those recently found by Arnold, Moore and Yaffe 
by analyzing the infrared divergences of the linearized collision integrals without screened 
interactions. Thus, this fact constitutes a non trivial check of the formalism we have used 
and also of the HTL approximation, because of the correct matching of the UV-divergences 
in our approach with the infrared divergences in the approach of Ref. |Q. 

With respect to the computation of the complete leading order of the transport coeffi- 
cients, we believe that, in analogy with the scalar theory, the introduction of all rungs made 
of the one- loop four-point functions may be of interest. Obviously, most of these rungs will 
lead to non divergent results, and only a few of them would yield the infrared divergences 
of the hard contribution matching with the logarithmic divergences we have found here. To 
carry out these computations within the framework we have presented, the classification of 
the four-point functions and the corresponding spectral densities would have to be studied. 
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APPENDIX A: THE SPECTRAL DENSITY OF THE RUNG IN A(/> 4 THEORY 



In this appendix, the spectral density in Eq. fl27|) is calculated. The product of two free 
Matsubara propagators is conveniently written as the double spectral representation 

duii du)2 1 



G (iu q + iiy k ,q + k)G (iiy k ,k) = J J 



(27r) 2 iv q - LU 2 
Po(u>i + w 2 , q + k)p (^i, k) p (uJi, q + k)p (^i - w 2 , k) 



(Al) 



iv k — U)\ iv q + iv k — U)\ 

where Po(p°,p) = 27r sgn(p°) <5((p ) 2 —p 2 )- By performing the sum over vj., and taking the 
imaginary part of the analytic continuation iv q —>■ q° + i0 + , one obtains 

d 3 h duj\ 



p(Q) 



po(u\ + q°, q + k)p (^i, k) - p (^i, q + k)p (^i - q°, k) 



(A2) 



The second piece of this integral may be rearranged by a shift of the uj\ integration, uj\ — > 
Ui + q°. Thus, after the use of the parity property p (cui+q°, q+k) = — p (— cj! — g°, — q— k)), 
and the substitution 

1 



p (-fc -g°,-q-k) 



Eq. 



yields the desired result 



p{Q) 



1 + n b (-k° - q°) 



G>(-K-Q), 



(27T) 



(A3) 



(A4) 



APPENDIX B: SUM RULES 



The required integrals over the Landau damping range of the frequency follow from the 
sum rules ]2l| derived from the analytic properties of the effective propagators. With the 



notation of Ref . |L6| , these sum rules are 

r« dq° 1 mi 



2vr q' 



™>d z L (q) 

q 2 (q 2 + m 2 D ) u L (q) 2 '' 



(Bl) 
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i dq° 



q°*PdQ) 



rn 



3g f-^ L (g), 



q 2ir 
i dg° 1 

-q^TT^ 



+ 



9q 2 
1_ _ z T (q) 
q 2 cj T (g) 2 ' 

1 - z T (q), 



ZL(q)uJL(q) 2 , 



q dq \ q y*pL(Q) 

l -o*MQ) 

£^(q°) 3 *MQ) = f + ^ - ^(gV T (g) 2 , 



(B2) 
(B3) 

(B4) 
(B5) 
(B6) 



where 2l,t(<z) are the residues at the quasi-particle poles. With the aim to extract the 
logarithmic dependence of the transport equations on q c , we use the the approximations 
valid at large q, q ^> thd, 

.2 / 2g 2\ 



Zh(q) 

Mqf 

u) T (q) 2 



m 



exp —2 — 



D 



rn 



1 + 



3m 2 



D 



Aq 2 



1 + - In 

3 



2 / ' 

D / 

1 . I ' m 2 



D 



1 + 4exp -2 



2g 



m 



m 



D 



+ 



m 



D 



Aq 2 



1 / m 
l + -ln 



8g 2 



(B7) 
(B8) 
(B9) 
(BIO) 



and we find the quadratic terms in m 2 D which give rise to the logarithmic terms ln(g c /m£>) 
after the q integration, 



rn 



D 



q 



4-2n 



X < 



n 4-2n 



X < 



1 n = 
1/3 n = 1 , 
1/5 n = 2 

-1/4 n = 
-3/4 n = 1 
-11/12 n = 2 



(Bll) 



(B12) 



The explicit form of the functions (3^ and /5t is 

Arrm 2 D q q9(q 2 - q 2 ) 



/5l(9o,c0 
/M^o, q) 



4g 2 {g 2 + m^[l - <2(g /g)]} 2 + tt 2 ^? 2 



Aq 2 {(g 2 - qf) [2q 2 + m 2 D Q(q /q)} + m 2 D q 2 } 2 + 7T 2 m 4 D q%(q 2 - gg) 



(B13) 
(B14) 
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where the function Q(x) is given by 



Q(x) = |lni±|. (B15) 



APPENDIX C: THERMAL WIDTHS 

To leading order in g (e) , the thermal width of a hard particle is obtained by the insertion 
of one HTL gluon (photon) propagator into the skeleton graph for the one-loop self-energy. 
For a hard quark (or a charged fermionic species s), the resulting expression is |16 

J4/ 



1? =^)^*-P'^ 



(CI) 



with g 2 C2(r) replaced by q 2 s e 2 in the case of a charged fermion. For a hard gluon, the same 
expression is valid if C^r) — > iV c . 

Next, we present the expressions for the thermal widths for hard particles to order 
g 4 lng -1 . The fermion self-energy is showed in Fig. ^, and is written as 

Z(tcu n ,p) = -g 2 C 2 (r)Tj2 / ^ 7 M *S(iu n + iu n , p + q) 7 ^(i^ q). (C2) 
Making use of the spectral representation of the soft fermion in the HTL approximation, 

^n,q) = / it- -7T— > (C3) 



oc 



27r g° — icj, 

with *A± verifying the parity properties 

*A ± (g ,q) = *A T (-g°,q), (C4) 

we may perform the Matsubara sum. After, one may extract the imaginary part of the 
continuation iu n — > p° + i0 + . When p° = p, the dominant contribution comes from the piece 
of the integrand which multiplies to 5{p° — q° — |p — k|). The replacement of this delta by 
6(q° — p ■ q), valid for g<p, selects the Landau damping piece of *A±, and the expansion 
of the remaining terms to lowest order in q/p yields 



1 

g 2 C 2 (r) (I \ r d 4 Q 



1? = ~ Tp tr Im S rct (p° = p, p)J = -IZ(p, ±)£ re V = p, p)u(p, ±) 



«*(p)) / 7^%°-p-q)^(Q)> (C5) 



2p \2 J J (2tt) 
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where 



T{Q) = (3 + {Q) 1 



+ 1 + 



(C6) 



The spectral functions (3±(Q) are given by 

^loliqT qo)0{q 2 - ql) 



^±(?o,q) 



7r2 ^pgo(gTgo) 2 ' 
4g 2 



(C7) 



[±<?o<?(<? =F ?o) + {±?o[i - Q(qo/q)] +qQ(qo/q)}] 2 

where the frequency plasma for the fermion is Uq = g 2 C2(r)T 2 /8 or q 2 e 2 T 2 /8. The inte- 
gral (p5|) has been treated in refs. 0, |(J and, for q c ^> u S) gives 



|«^-p. qOT Q) = g 



hi 



V2< 



l + ln2 



u s 



(C8) 



The 0(<7 4 lng _1 ) thermal width for a hard gauge boson is associated with the imaginary 
part of the self-energy of the diagram in Fig. [?]. It reads 

d 3 q 

where the prefactor 2 comes from the two possible ways to arrange a soft fermion propagator 
in the graph. A similar treatment to the fermionic case leads to the result 



/c^q 
— — - tr {7^ S(iu n + iu n , p + q)j u *S(iis n , q)} , (C9) 



= 9 ^f^ (~ ~ n f <pj) J «6tf - P • *)HQ), (CIO) 

where IIt(P) = p l p>)/2. For the case of a hard photon, the correct result is obtained 

with the substitutions g 2 NfC(r) — > e 2 ^charged species Q 2 an d —* ^s- 
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